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Abstract 

In this paper we analyze Yukawa unification in a three family SO(IO) SUSY 
GUT. We perform a global analysis and show that SUSY effects do not de¬ 
couple even though the universal scalar mass parameter at the GUT scale, rriiQ, 
is found to he between 15 and 30 TeV with the best fit given for rriiQ ~ 25 TeV. 
Note, SUSY effects don’t decouple since stops and bottoms have mass of order 5 
TeV, due to RG running from Mqut- The model has many testable predictions. 
Gauginos are the lightest sparticles and the light Higgs boson is very much Stan¬ 
dard Model-like. The model is consistent with flavor and GP observables with 
the BR{y, —)■ ey) close to the experimental upper bound. With such a large 
value of mi6 we clearly cannot be considered “natural” SUSY nor are we “Split” 
SUSY. We are thus in the region in between or “SUSY on the Edge.” 
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1 Introduction 


Gauge and Yukawa unified SO(IO) supersymmetric grand unified theories (SUSY GUTs) 
m are a class of highly constrained models. By performing a global £t to low energy data, 
such as fermion masses and mixing angles and flavor observables, SO(IO) SUSY GUT models 
can provide experimental constraints on sparticle masses accessible at the LHG energy scale. 
The specihc SO(IO) SUSY GUT model that we study in this paper has a [U(l) xZ 3 xZ 3 ] 
family symmetry [HIE!. This model includes three families of quarks and leptons and has 
been shown to provide good hts to low energy precision electroweak observables, including 
fermion masses and mixing angles [HHI2]. 

In this paper we extend the analysis of previous papers by including additional low 
energy flavor observables. We also identify which observables are better fit by incorporating 
the SUSY effects. We then discuss the many experimental tests of this model. Finally we 
evaluate the amount of hne-tuning in this scenario. As might be expected the amount of 
fine-tuning is enormous. However, we show that introducing certain boundary conditions at 
the GUT scale has the effect of dramatically reducing the amount of hne-tuning. The origin 
of such boundary conditions in this theoretical sweet spot is not known, although we discuss 
some possible sources in an Appendix. 

The paper is organized as follows. In Section]^ we present an overview of the SO(IO) 
model and matorQ, the program that runs all the parameters and calculates most of the 
observables. In Section we present the results of a global analysis for both the universal 
and “mirage mediated” gaugino mass boundary conditions. We show that is minimized for 
values of the universal scalar mass, rriiQ of order 25 TeV, for both universal gaugino masses 
at the GUT scale or “mirage mediation” with non-universal gaugino masses, respectively 
(see Figs. and |^. The function now includes the GP-even angular observables of 
measured by LHGb. Some of these observables are in slight tension with the 
standard model (SM) prediction [T31 [Hj. In particular, the observable Pg integrated over 
the dimuon invariant mass squared, 1 < < 6 GeV^, has a 2.5a discrepancy from the SM 

prediction. In addition, the angular observable P 4 integrated over 14.18 < < 16 GeV^ 

also has a 2.7a tension with the SM prediction. We also compare our hts to inclusive vs. 
exclusive measurements of the GKM elements Wb and Wb- 

The hrst two families of squarks and sleptons have mass of order rriiQ and thus they 
decouple from low energy physics. However, the third family scalars are significantly lighter, 
with mass between 4 and 10 TeV. This is a consequence of a natural inverted scalar mass 
hierarchy [15]. Hence they do not decouple and thus contribute to lowering In Section 3.2 


we discuss the observables which are better ht by including SUSY loops. We then discuss 


bounds on the gluino mass in Section 3.3 and further predictions in Section 3.4 


Since our scalar masses are heavy there is, in general, considerable hne-tuning of order 
1 part in 10^. In Section we present the results of a hne tuning analysis of the model 
and demonstrate that given particular ratios of soft breaking parameters at the GUT scale 
the amount of hne tuning is minimized to roughly 1 part in 500 (perhaps these boundaries 
conditions can be obtained in a fundamental SUSY breaking model). Possible origins for the 


^ “maton” was developed by Radovan Dermisek and modified for large scalar masses by Archana Anan- 
drakrishnan, B. Charles Bryant, Zijie Poh and Akin Wingerter. 
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necessary boundary conditions are discussed in the Appendix. 

Finally, we note that our soft breaking parameter rriie lies in the range between 15 and 
30 TeV. Hence the gravitino mass is expected to be large, perhaps large enough to avoid 
the cosmological gravitino problem [16]. The moduli may also be suitably heavy to avoid a 
cosmological moduli problem [T7M20] . At the same time we see that SUSY effects have not 
decoupled from low energy flavor observables, thus our soft SUSY breaking parameters lie 
on the edge between Split and natural SUSY. We denote this region of parameter space as 
SUSY on the Edge. Gauginos are expected to be visible at the LHC with gluinos decaying 
predominantly into third generation quarks with special decay signatures. While the light 
Higgs discovered by CMS and ATLAS is expected to be very much Standard Model-like. We 
also present predictions for additional flavor and CP violating observables, see Table In 
Section]^ we conclude and discuss phenomenological tests of our model. 


2 Model and Procedure 


The details of our SO(IO) SUSY GUT model is presented in [9|. The GUT scale boundary 
conditions that we use are universal squark and slepton masses, mie, universal trilinear 
couplings, Aq, and non-universal Higgs masses, and As for the boundary condition 
for the gaugino masses, we consider two different cases. One with universal gaugino masses, 
Mi/ 2 , and another with mirage mediated gaugino masses. The mirage mediated gaugino 
masses are dehned as [10] 





g^b^a MpA 

miQ J 


Mi/2 , 


( 1 ) 


where M 1/2 and a are free parameters and bi = (33/5,1, —3) for z = 1, 2, 3 are the appropriate 
/3-function coefficients. Mirage boundary condition is interesting because a = 1.5 is the 
optimal scenario for a well-tempered dark matter [21] . On the other hand, it has been shown 
that the LSP of universal boundary condition is predominantly Bino-like m, which leads 
to an over-closed universe. However, this problem can be solved by introducing axions with 
mass lighter than the LSP into the model. Hence, in this paper, we study the cases where 
the universal boundary condition, a = 0, and mirage boundary condition with a = 1.5. 

On the other hand, the charged-fermion sector of our model has 12 parameters - 11 
Yukawa parameters and tan/3. In comparison, the SM has 13 free parameters. Hence, our 
model has 1 prediction in the charged fermion sector. Including the neutrino sector, our 
model has 3 additional free parameters to £t 6 observables - 2 Am^s, 3 real mixing angles 
and one CP violating phase. Hence, our model has 4 predictions in the fermion sector. 

To summarize, all the free parameters of our model are listed in Table With universal 
gaugino masses boundary condition, our model has 24 free parameters while with mirage 
mediated gaugino mass boundary conditions the number of free parameters increases to 25. 

The detailed procedure of our calculation is presented in [9]. In addition, we want to 
emphasize that after integrating out the right-handed neutrinos, we use the two-loop Minimal 
Supersymmetric Standard Model (MSSM) renormalization group equation (RGE) to run 
down to the weak scale. At the weak scale, we include the one-loop threshold corrections 
for the Yukawa and the gauge couplings, which was calculated by Pierce et. ah [22]. The 
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Sector 

Universal Gaugino Masses 

No. 

Mirage Mediated Gaugino Masses 

No. 

Gauge 

UG, Mq, £3 

3 

no. Mg, £3 

3 

SUSY (GUT scale) 

mi6, Mi/2, Ho, uih^, 

5 

ruie, Mi/2, Hq, mn^, mn^, a 

6 

Yukawa Textures 

e, 6 , A, p, fJ, e, (pa--, (pei 

11 

P, fT, £, (j)p, pa, Pi, Pi 

11 

Neutrino 

Mr,, Mr^, Mr,, 

3 

Mr,, Mr^, Mr„ 

3 

SUSY (EW Scale) 

tan P, /i 

2 

tan/ 3 , M 

2 

Total 


24 


25 


Table 1: Our model has 24 free parameters for universal gaugino masses boundary condi¬ 
tions. For mirage mediated gaugino masses boundary condition, there is an additional input 
parameter, a, which determines the amount of splitting of the gaugino masses at GUT scale. 


gluino mass. Mg, and CP odd Higgs mass, Ma, are pole masses. However, we do not include 
one-loop threshold corrections for the other scalar masses. Instead, we estimated that these 
corrections are about 10%. Hence, the GUT Scale parameters of our model have an inherent 
10% theoretical error. Adding the one-loop threshold correction for the soft scalar masses 
can be a future project. 

In addition to calculating the observables included in [9j, we included a low bin (1 < 

< 6 GeV^) and a high g^ bin (14.18 < g^ < 16 GeV^) for each of the following 4 
CP-even B —)■ angular observables: Fi^,P 2 ,Pl, and Pg. We did not include other 

CP-even observables because the theoretical uncertainty of those observables are much too 
big. Hence, they do not constrain our model. These angular observables are calculated by 
superiso version 3.4 [23]. Since superiso assumes that all soft parameters are real and only 
takes the diagonal entries of the trilinear couplings into account, we do not include CP-odd 
B —)■ A*/i/i angular observables in our analysis. 

Table includes 45 observables that we include in our global analysis. The pro¬ 
gram that calculates these observables and the theoretical errors that are assigned to each 
observable are also included in the table. The theoretical uncertainty for the B —)■ ob¬ 

servables are taken from the superiso manual. However, since superiso does not take into 
account the imaginary part of the soft parameters, we assigned an additional 15% theoretical 
errors to the calculation. 

To perform the global analysis, we construct a x^ function 


X 


2 


E 


I ^th 
\^i 


- X 


exp 12 


a" 


( 2 ) 


where xf' are the calculated values, are the experimentally measured values, and are 
the sum of the squares of the experimental and theoretical uncertainties, which are also listed 
in Table To hnd the minimum of this x^ function, we use the Minuit package maintained 
by CERN [23]. As in most minimization problems, obtaining the true global minimum is 
not guaranteed. To increase the likelihood of obtaining the true global minimum, we iterate 
the minimization process with random initial guesses for the free parameters. 

To calculate the x^/d.o.f, we assume that the observables are uncorrelated. So, for the 
universal boundary condition, we have 45 — 24 = 21 degrees of freedom, while for the mirage 
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boundary condition, we have 45 — 25 = 20 degrees of freedom. Given these gross assumptions, 
one should not take the value of x^/d.o.f too seriously. 

3 Results: Global ^ Analysis 

Several benchmark points with the results of the global analysis are given in Appendix 
[A| and a plot of as a function of the parameter mig is given in Figure [I] or contours 

in the two dimensional plane of Mg vs. irtiQ is given in Figure Let us now discuss some 

features of the results. 

3.1 Inclusive vs. Exclusive \Vub\ and \Vcb\ 

Due to the discrepancy between the values of |14b| and \Vch\ determined from inclusive and 
exclusive semi-leptonic decay. We dehne three different functions: 

1 . \Vub\ and \Vcb\ are taken to be the inclusive values 

2 . \Vub\ and \Vcb\ are taken to be the exclusive values 

3. \Vub\ and |Wb| are taken to be the average of inclusive and exclusive values with error 
bars overlapping with the error bars of both the inclusive and the exclusive measure¬ 
ments 

The results of these three analyses are shown in Figure [T] We see that for both the 
universal boundary condition a = 0 and mirage boundary condition with a = 1.5, the 
X^/d.o.f obtained by htting to the inclusive values are the biggest. Hence, we predict that 
the exclusive values of | 14 fe| and |Wb| are the correct values for both universal and mirage 
gaugino masses. 

Since the difference between case ( 2 ) and case (3) is small and to be conservative, the 
analyses of the rest of the paper are done for case (3), where |14fe| and |14b| are the average 
of the inclusive and exclusive values. 

3.2 SUSY Non-decoupled observables 

B physics observables 

Some of the measured angular observables of i? —)■ K *are in tension with the 
SM prediction. For example, P'^ in the high bin (14.18 < < 16 GeV^) has a 2.7a 

discrepancy with the SM prediction, in the low bin (1 < < 6 GeV^) has a 2.5cr 

discrepancy with the SM prediction, and P 2 in the low q^ bin has a 2a discrepancy with the 
SM prediction usi m EH. These observables are dehned in |32lE3]- la addition, previous 
analysis [aHEl] found that the tension in P 4 of the high q^ bin cannot be explained by the 
MSSM. On the other hand, the tension of Fl and P 5 of the low bin can be explained by 
the MSSM by having a negative contribution to the Gy Wilson coefficient. In the standard 
model Gy ~ —0.32. The tension in Fl and P5 can be further reduced by making Gy more 

negative [311 EH]- 
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Observable 

Exp. Value 

Ref. 

Program 

Th. Error 

Mz 


91.1876 ± 0.0021 GeV 

[25] 

Input 

0 .0% 

Mw 


80.385 ± 0.015 GeV 

[25] 

mat on 

0.5% 

Ckem 


1/137.035999074(44) 

m 

mat on 

0.5% 

G, 


1.1663787(6) x 10“® GeV“^ 

[25] 

mat on 

1 % 

oi-AMz) 


0.1185 ± 0.0006 

[25] 

mat on 

0.5% 

Mt 


173.21 ± 0.51 ± 0.71 GeV 

[25] 

mat on 

0.5% 



4.18 ±0.03 GeV 

[25] 

mat on 

0.5% 

Mr 


1776.82 ±0.16 MeV 

[25] 

mat on 

0.5% 

mb — me 


3.45 ± 0.05 GeV 

[25] 

mat on 

10 % 

mc{mc) 


1.275 ±0.025 GoV 

[25] 

mat on 

0.5% 

m,(2 GeV) 


95 ± 5 MeV 

m 

mat on 

0.5% 

msjmd (2 GeV) 


17-22 

m 

mat on 

0.5% 

Q 


21 - 25 

[25] 

mat on 

5% 

M^ 


105.6583715(35) MeV 

[25] 

mat on 

0.5% 

Me 


0.510998928(11) MeV 

[25] 

mat on 

0.5% 

IK.I 


0.2253 ± 0.0008 

[25] 

mat on 

0.5% 

|Vc6| (Inclusive) 


0.0422 ±0.0007 

[25] 

mat on 

0.5% 

|Vc6| (Exclusive) 


0.0395 ± 0.0008 

m 

mat on 

0.5% 

|V,i| (Both) 


0.0408 ±0.0021 

m 

mat on 

0.5% 

IKil (Inclusive) 


0.00441 ± 0.00024 

[25] 

mat on 

0.5% 

\Vub\ (Exclusive) 


0.00328 ± 0.00029 

[25] 

mat on 

0.5% 

|V„ 5 | (Both) 


0.00385 ± 0.00086 

[25] 

mat on 

0.5% 

\VtA 


0.00840 ± 0.0006 

[25] 

mat on 

0.5% 

\Vu\ 


0.0400 ± 0.0027 

[25] 

mat on 

0.5% 

sin 2/3 


0.682 ±0.019 

m 

mat on 

0.5% 



(2.2325 ± 0.0155) x 10“^ 

[25] 

susvf lavor[2^ 

10 % 

AmBjAniBj 


35.0345 ± 0.3884 

[25] 

susvf lavor[2^ 

20 % 



(3.337 ±0.033) x 10““ MeV 

[25] 

susvf lavor[2^ 

20 % 

Am2i 


(7.02 - 8.09) X 10-5 eV^ (3 ct range) 

[27] 

mat on 

0.5% 

Am§i 


(2.317 - 2.607) X 10“^ eV^ (3cr range) 

[27] 

mat on 

0.5% 

sin^ 9\2 


0.270 — 0.344 (3(7 range) 

[27] 

mat on 

0.5% 

sin^ 023 


0.382 — 0.643 (3(T range) 

1271 

mat on 

0.5% 

sin^ 9 i3 


0.0186 - 0.0250 (3t7 range) 

m 

mat on 

0.5% 

Mh 

125.7 ±0.4 GeV 

m 

split suspect |28^ 

3 GeV 

BR(& —> 57 ) 


(343 ± 21 ± 7) X 10“® 

m 

superiso|23l 

40% 

BR{Bs —> 


(2.87;/) X 10“® 

d] 

susvf lavor[2^ 

20 % 

BR(_Bd 


(3.9t!;h X 10“'“ 

[lo] 

susvf lavor[2^ 

20 % 

BR{B —> Ti/) 


(114 ±22) X 10-® 

[29] 

susvf lavor[2^ 

50% 

BR(B ^ K'ti+ti- 

)l<q2<6 GeV2 

0.34 ± 0.03 ± 0.04 ± 0.027:!55 x 10“^ 

[13 

superiso|23l 

105% 

BR(B ^ K'ti+ti- 

)l4.1S<g2<i6 GeV^ 

0.45 ± 0.06 ± 0.04 ± 0.047:“ x 10“^ 

El 

superiso|23l 

190% 

«ohAFB(B ^ K'fi 

V-)) 

4.9 ± 0.9GeV^ 

El 

superiso|23l 

25% 

Fl(B ^ K'ti+n- 

)l<p<6 GeV^ 

ri c'r+0.08+0.03 

07-0.03 

E3 

superisol23l 

45% 

Fl{B K* fl )i4.1S<g2<i6 GeV2 
— 2 P 2 — A^‘^{B —> Qey2 

/-V nQ+0.08+0.03 

U.00_o.o7-0.03 
n ec+0.24+0.04 

U.DO_o, 22 - 0.01 

E3 

E3 

superiso|23l 

superiso|23l 

80% 

95% 

- 2 P 2 = A^‘^{B P*M''’M”)l4.18<g2<i6 GgV 2 

1 nn+000+0.00 

i-'JU_o.o5-0.02 

E3 

superiso|23l 

45% 

PUB ^ 

l<g2<6 GeV2 

0.58172 ± 0.06 

E3 

superiso|23l 

30% 

Pl(B ^ 

14.18<q2<16 GeV^ 

-0.187:S±0.08 

E3 

superiso|23l 

35% 

Pi{B K’ti+n- 

l<g2<6 GeV^ 

o.2i7;“ ± 0.03 

Ea 

superisol23l 

45% 


14.18<q2<16 GeV^ 

-0.797:i3±0.18 

Ea 

superiso|23^ 

60% 


Table 2: This table includes 45 observables that we fit. All experimental errors are la unless 
otherwise indicated. Column 4 shows the software package that gives us the theoretical prediction. 
Mz is fit precisely to impose electroweak symmetry breaking. To account for the inconsistencies 
in the inclusive and exclusive measurements of \Vub\ and \Vcb\-, we perform the global analysis 
using the inclusive and exclusive measurements separately. We also perform an additional analysis 
where the error bars of \ Vub\ and 11461 cover both the inclusive and exclusive values. In additional 
to the theoretical errors indicated in superiso manual, we added an additional 15% error to each 
of the B —)• K *observables because superiso does not take into account the phases of soft 
terms. 
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Figure 1: This plot shows the value of x^/d.o.f as a function of rriiQ for cases where the value 
of 1 14 b I and 1 14 b I are taken to be the inclusive values, the exclusive values, or the average of 
inclusive and exclusive values. Solid lines refers to the universal boundary condition, a = 0, 
while dashed lines refer to the mirage boundary condition with a = 1.5. This plot shows 
that our model favors the exclusive values of |14b| and |14b|- 


In the MSSM, chargino-stop loops and charged Higgs loop contribute to Cy. The Cj 
contribution from the charged Higgs is always negative. The charged Higgs of our model has 
mass around 2 TeV. So, the charged Higgs contribution to Cj is non-negligible and is in the 
correct direction. 

The chargino-stop loop contribution of has the following form [Sj^ 

mssm ^ M,tan4 ^.g^ _ ^3^ 

ml 

Since sign(/iH4 is negative in our model, this term contributes to Cj in the wrong direction. 
Hence, to reduce the contribution of this term, our model favors large scalar masses. 

From our global analysis, we see that the calculated value of P'^ in the high bin does 
not depend on myg, which is expected. In addition, the value of P'^ calculated in our model 
is in agreement with the SM. Hence, our results are in agreement with previous analysis that 

^Note, the equation for the chargino contribution to gjygjj Eqn. 21, Ref. [3T] apparently has 

the wrong sign. 
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the tension in cannot be explained in the MSSM. As shown in Table the tension of Fl 
and P 5 with the experimental values decreases as mig increases. This is again in agreement 
with our expectation as explained above. 

SUSY corrections to the W mass 

On the other hand, the correction for is given by |36f[38] 


(5M, 


M^ 


w 


w 


-w 


r2 _ s2 


Ap 


(4) 


and the 1 -loop squark contribution is given by 
3G,. 




8\/2: 


TT^ 

.2„2 ; 


l-s^-c^-Fo(m^^,mJ - ,mj 


bi’ 


62 ^ 


(5) 


+ c^CrFoimr ,m-J + c^SrFo{mr , m-J + s^cpFo{m^ m- ) + s,-SrPo(mr , 


( 6 ) 


where sw = sin 6 'iy, cw = cosOw, Sq = sin 6 ^ 5 , Cq = cosOq, and 

17 f \ ^ 2xy X 

Po(x, y) = x + y -In - . 

x-y y 


(7) 


pQ has properties of Po(x,x) = 0 and Po(x, 0) = x. Hence, we see that when the mass 
splitting of the squarks is large, the SUSY contribution to the 1-loop Mw can be signihcant. 
This is in agreement with our analysis which shows that the pull from Mw increases as the 
value of mi 6 increases above 20 TeV. Hence, SUSY corrections to are signihcant and 
they can go in the right direction. 


Light Higgs mass 

Fitting to the Higgs mass also constrains the value of miQ. The dominant one-loop 
contribution to the Higgs mass is given by 


ml ~ 


m| cos^ 2/3-1- 


m; 


(47r)2 


/\/f2 


+ 


X 


mt 


^''^SUSY 


X? 


12MgijgY 


( 8 ) 


where Xt = At — pj tan /3 is the stop mixing parameter and Mg^gY = model, 

Xt < —\/6 Msusy and the ratio Xt/MsusY becomes less negative as mig increases. Hence, 
as mi 6 increases, the Higgs mass also increases. The pull in due to has a minimum 
around mie = 25 TeV. 

Hence, the contributions of Mw, Mh, and 6 -physics observables to as listed in Table 
1^ help explain the shape of as a function of (see Figure [^. 












Pull 

mi6 

10 

15 

20 

25 

30 

M\y 

0.2110 

0.1878 

0.1851 

0.2320 

0.3981 

Mh 

2.5474 

1.1795 

0.3454 

0.1882 

0.6582 

BR(B — )■ Tu) 

1.1978 

1.3952 

1.3557 

1.3588 

1.3771 

Fl{B A'*/i+/i“)i<g2<6 GeV2 

0.2696 

0.2488 

0.2219 

0.2101 

0.2057 

Pi{B —)■ A'*p+/i“)i<q2<6 GeV^ 

1.7066 

1.7066 

1.7066 

1.7066 

1.7066 

t A'*/i+/i“)i<Q2<6 GeV^ 

2.4110 

2.3432 

2.2746 

2.2451 

2.2339 


14.1511 

8.9744 

7.2154 

7.0206 

7.5220 


mi6 

Fit Value 

Exp. 

Value 

10 

15 

20 

25 

30 

Mw 

80.4699 

80.4606 

80.4595 

80.4784 

80.5454 

80.3850 

Mh 

117.9901 

122.1303 

124.6547 

126.2697 

127.6920 

125.7000 

BR{B tv) X 10 ^ 

6.6329 

6.1340 

6.2299 

6.2223 

6.1778 

11.4000 

Fl{B — )■ K* GeV^ 

0.7434 

0.7353 

0.7251 

0.7207 

0.7191 

0.6500 

Pi{B — )• ^ )l<q^<6 GeV^ 

0.8174 

0.6711 

0.5921 

0.5717 

0.5657 

0.5800 

Pi{B — )• K*^'^^ )l4.18<g2<i6 GeV^ 

1.2190 

1.2190 

1.2190 

1.2190 

1.2190 

- 0.1800 

P^iB —>■ K*^^'^M)l<q^<6 GeV^ 

- 0.7301 

- 0.5529 

- 0.4625 

- 0.4335 

- 0.4235 

0.2100 


Table 3: This table shows the set of observables with a = 0 and Mg 1.2 TeV. As argued 
by Altmannshofer et. al. B —)■ favors large [3T], which is in agreement with our 

analysis. On the other hand, as explained in the text, htting to and disfavors large 
m^. These effects collectively contribute to having a minimum around mig = 25 TeV. 
The value is plotted as a function of mie in Figure In addition to the observables 
contributing directly to a minimum of at = 25 TeV, also included in this table is 
the calculated value of at high bin. This is to illustrate that the value of does not 
depend on the mie, which again agrees with previous analysis [HI IM] . 
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mi6(TeV) 


Figure 2: This figure shows the contribution to as a function of mig just from the set of 
observables listed in Table Fitting to the values of Mw, Mh, and 6-physics observables 
listed in Table helps explain why is minimized at mig ~ 25 TeV. 


3.3 Bounds on Mg 

To obtain a better picture for the favored value of gluino mass, we plotted two contour plots 
of Mg vs. rriiQ. One for the universal boundary condition a = 0 and another for the mirage 
boundary condition with a = 1.5. The current bound on Mg, for our model, is around 
1.2 TeV [12]. The contour plots are created by calculating y^/d.o.f for 25 equally distributed 
values of mig and M 1 / 2 , which gives us 1 < Mg < 2 TeV. We then use cubic interpolation 
to obtain the smooth contours of x^/d.o.f. 

In addition to the contour lines of the y^/d.o.f, we also plotted a 4 ct contour line. From 
this, we see that for mirage boundary conditions Mg < 1.8 TeV. However, for universal 
boundary condition, the Aa Mg bound can be as high as 3 TeV, which is not shown in 
the Figure. Hence, for mirage boundary conditions, we expect the 4cr bound on the gluino 
mass to be within reach in the next run of the LHC. As pointed out by na, the dominant 
decay mode of the gluino in the universal gaugino mass boundary condition is t b xf. The 
remaining decay modes are t t and 6 6 x° for i = 1,2. On the other hand, the dominant 
mode for gluino decay in the mirage gaugino mass boundary condition is t b xf for i = 1,2. 
In all cases, the dominant signature for gluinos in this model is given by 6 jets, leptons and 
missing Et m 
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X^/d.o.f Contour Plot with a = 0 



mislTeV) 


. TTl 

2.50 2.75 3.00 3.25 3.50 3.75 4.00 4.25 



r .Ill 

2.4 2.7 3.0 3.3 3.6 3.9 4.2 4.5 4.8 

Figure 3: These plots show the contour of x^/d.o.f as a function of the Mg and rriiQ. The 
4cr bound is also included in the plot. For a = 1.5, the upper bound is within reach of the 
next run of LHC. In addition, we also see that our model favors ruig ~ 25 TeV. 
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3.4 Additional predictions for some benchmark points. 

The mass spectrum of the benchmark point of Mg 1.2 TeV and mie = 25 TeV is shown 
in Table From the analysis, we see that the scalar masses are predicted to be around 
5 TeV, while the hrst and second generation scalars have mass around mig ~ 25 TeV. With 
scalars in this mass range, the stop in our model does not completely decouple and can have 
non-negligible effects on flavor physics. In addition, our light Higgs is SM-like with the heavy 
Higgs with mass around 2 TeV. 

In Table we give the light sparticle masses, the CP violating angle for neutrino oscilla¬ 
tions, 5, the branching ratio BR{^ —)■ ey) and the electric dipole moment of the electron for 
two different values of Mg and for a = Oandl.5. Note that, in general, the gauginos are the 
lightest sparticles. In addition, BR(/j, —)■ ey) and the electric dipole moment of the electron 
are within reach of future experiments. 


mi6 

25 

25 

25 

25 

a 

0 

1.5 

0 

1.5 

xVd.o.f 

2.158 

2.275 

2.220 

2.505 


4.903 

5.011 

4.909 

5.249 


6.021 

6.120 

6.033 

6.301 

rn-b. 

5.989 

6.088 

6.455 

6.606 


6.454 

6.541 

6.445 

6.267 


9.880 

9.931 

9.912 

10.040 

rrif^ 

15.369 

15.365 

15.393 

15.516 

Mg 

1.202 

1.187 

1.613 

1.690 

m^o 

0.203 

0.551 

0.279 

0.900 


0.404 

0.665 

0.538 

1.018 

^xt 

0.404 

0.665 

0.538 

1.018 

^xt 

1.128 

1.243 

1.232 

1.537 

Ma 

2.194 

2.082 

2.477 

3.352 

sin (5 

-0.289 

-0.482 

-0.520 

-0.576 

BR{^ —)■ ey) x 10^^ 

1.108 

1.430 

1.239 

1.340 

edmg X 10^° (e cm) 

-1.403 

-3.305 

-1.763 

-5.886 


Table 4: Predictions with rriie = 25 TeV for Mg 1.2 TeV and 1.6 TeV. All masses in the 
table are in TeV units. Our prediction for the branching ratio /r —>■ ey is consistent with the 
current upper bound of 5.7 x 10“^^ [25]. In addition, our prediction of the electron electric 
dipole moment is consistent with the current upper bound of 10.5 x 10“^® e cm [23] . 
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4 Results: Fine-Tuning 


We studied the fine-tuning of our model using the fine-tuning measure introduced by Ellis 
et. ah [39], and studied in detail by Barbieri and Giudice [lO] . 


Abg = max A^,, A^^ 


a In M| 
dlna^ 


(9) 


where OjS are input parameters of the model. This fine-tuning measure calculates the sensi¬ 
tivity of Mz due to a small variation of the input parameters defined at GUT scale. 

Electroweak symmetry is broken radiatively in our model. From radiative electroweak 
symmetry breaking, the GP-odd Higgs mass, tua, and the /i-term are calculated at one-loop 
[22] • This calculation requires the physical Z pole mass, Mz- Hence, in our model, Mz is 
fit precisely. To make sure that radiative electroweak symmetry breaking is consistent, tua 
and fi are calculated iteratively until they converge. 

On the other hand, when we calculate fine-tuning using ([^, we use the benchmark points. 
The benchmark points are the inputs that produce minimum value for their respective 
value of rriiQ and Mi/ 2 - Hence, at each benchmark point, radiative electroweak symmetry 
breaking is consistent. Thus, instead of fixing Mz and calculating rriA and p iteratively, we 
then use the value of niA and /i to calculate Mz- We then compare this value of Mz to the 
exact value to obtain the fine tuning parameter Abg- 

The input parameters that we vary are a* = {/i, M 1 / 2 , mie, Hq}. The results of 

our calculations are summarized in Table These results can be understood by the running 
of GUT scale parameters that contribute to Z mass. For tan/3 = 10, Mz written in terms 
of GUT scale parameters is [HHii] 


M| -2.18p2 + 4:-22M^/2 - O.82M1/2H0 + 0.22^2 - 1.27m^^ - 0.053m^^ + 1.34m?6 • 

( 10 ) 


Although the above equation is derived for tan (3 = 10, to the lowest order approximation, we 
do not expect this result to change drastically when tan (3 increases to ~ 50. The calculated 
fine tuning values shown in Table are of the same order as the fine tuning predicted from 
this equation. As an example, by direct substitution of the rriie = 20 TeV benchmark points 
into the above equation, we find that if and Ao/mie are fixed, then Abg ~ 2000 

is of the same order as our calculation. 

From Table we see that if there are no constraints on the input parameters (first five 
rows), then the fine-tuning is about 1 part in 10 ^. However, if the GUT scale parameters 
are constrained such that ~ y/2 and Ao/miQ ^ —2, then the fine-tuning of our 

theory is reduced to about 1 part in 500 . This suggests that, in a more fundamental natural 
theory, the ratio of rriiQ with ^ and Aq could be fixed naturally. In the context of the Bear 
et. ah argument that one should combine dependent terms into a single independent quantity 
before evaluating fine-tuning|TT], we claim that and Aq might be dependent 

quantities in a more fundamental theory. Hence, one should combine these quantities before 
calculating fine-tuning. We discuss one possible partial example in the Appendix. 
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Fine-Tuning of Benchmark Points with a = 0 and Mg 1.2 TeV 


Varying Parameters 

mi6 

lOTeV 

15TeV 

20TeV 

25TeV 

30TeV 

fi 

140 

190 

210 

360 

490 

Mi/2 

260 

340 

400 

430 

450 

mi6 

12000 

27000 

47000 

74000 

110000 


760 

1500 

3900 

6100 

8700 


10000 

23000 

40000 

62000 

89000 

^0 

9300 

21000 

39000 

61000 

85000 

mi6 with Ao/rriiQ hxed 

22000 

49000 

87000 

130000 

190000 

mi6 with niH^^d/mie hxed 

9500 

22000 

40000 

62000 

86000 

mi6 with mH^^/miQ,Ao/mie hxed 

240 

400 

630 

740 

850 


Fine-Tuning of Benchmark Points with a = 1.5 and Mg 1.2 TeV 


Varying Parameters 

mi6 

lOTeV 

15TeV 

20TeV 

25TeV 

30TeV 

fi 

110 

240 

290 

340 

380 

Mi/2 

320 

420 

500 

560 

580 

mi6 

12000 

28000 

48000 

74000 

110000 


750 

1500 

4600 

6000 

8500 

rriHu 

10000 

23000 

39000 

62000 

89000 

^0 

9200 

21000 

39000 

60000 

86000 

rriiQ with Aq/itiiq hxed 

22000 

49000 

87000 

130000 

190000 

mi6 with hxed 

9600 

21000 

39000 

61000 

87000 

mi6 with mH^ ,^/mi(i,Ao/mi(i hxed 

330 

450 

670 

890 

1100 


Table 5: Without hxing any ratios, the hne-tuning is 1 part in 10®. When the ratio of 
niH^a/'miQ and Ao/mig are hxed, the hne-tuning according to (|^ is about 1 part in 500. 
This suggests that these ratios should be hxed in a more fundamental natural theory. In 
addition, hne-tuning increases as mig increases. Hence, in our model, small rriiQ is favored 
by naturalness. 
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5 Conclusion: SUSY on the Edge 


We have analyzed a three family SO(IO) SUSY GUT with Yukawa unihcation for the third 
family. The model gives reasonable hts to fermion masses and mixing angles, as well as many 
other low energy observables; see Appendix with some benchmark points of the global 
analysis. A plot of as a function of the parameter rriiQ is given in Figure [I] or contours 
in the two dimensional plane of Mg vs. mie is given in Figure]^ 

We performed an analysis with universal gaugino masses and with non-universal gaugino 
mass with splitting determined by “mirage mediation” boundary conditions described in 
Eqn. The parameter a = 0 for universal gaugino masses and we also take a = 1.5 which 
is consistent with a well-tempered dark matter candidate [2T] . In both cases the model favors 
mi6 ~ 25 TeV. Nevertheless, due to RG running na, stops and sbottoms have mass of order 
5 TeV, while the hrst two family scalar masses are of order mig. With tjiiq lying in this mass 
range, stops in our model do not completely decouple from low energy flavor observables 


(see Section 3.2). Best hts are found with a gluino mass less than 2 TeV. Our gluinos decay 

Moreover, in a previous analysis 


predominantly into third generation quarks [I^. Moreover, m a previous analysis [i^ we 
showed that the dominant LHG signature for gluinos in the model is given by 5-Jets, leptons 
and missing Et- Note that, in general, the gauginos are the lightest sparticles. The GP 
odd Higgs mass is of order 2 TeV, thus the light Higg couplings are very much Standard 
Model-like. In Table we present additional predictions. We give the predictions for the 
GP violating angle for neutrino oscillations, 6, the branching ratio BR{fi —)■ ey) and the 
electric dipole moment of the electron for two different values of Mg and for a = Oandl.5. 
In addition, BR{^ —)■ ey) and the electric dipole moment of the electron are within reach of 
future experiments. Thus this theory is eminently testable! 

We evaluated the amount of high scale hue-tuning of our model. In general we hnd hne- 
tuning of order 1 part in 10^. However we note that with particular boundary conditions 
at the GUT scale (when the ratio of mig to Aq and are hxed at Ao/'^^ie ~ ~2 and 

^/mi6 ~ \/2) the hne-tuning is reduced to 1 part in 500. We do not have a fundamental 
theory that gives these two ratios naturally. Nevertheless, in such a fundamental theory the 
amount of hne-tuning is reduced considerably. 

Finally, with the large value of mie ~ 25 TeV we expect the gravitino mass to be at 
least this large. Perhaps it is large enough to avoid a cosmological gravitino problem [16]. In 
addition, moduli may also be suitably heavy to avoid a cosmological moduli problem HMD]. 
Hence the scalar masses are clearly in an intermediate range, i.e. too heavy to be “natural” 
and lighter than “Split SUSY.” We thus are positioned on the border between these two 
limiting cases, i.e. this is “SUSY on the Edge.” 
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A 


Benchmark Points 


Table 6: Benchmark point with mie = 25 TeV, Mg = 1.202 TeV, a = 0: 
(l/acMcea) = (25.98,2.55 x IQi® GeV,-1.30%) 

(A, Ae, cr, Ae, p, Ae', Ae^ = (0.6101,0.0308,1.1559,0.0049,0.0698, -0.0019,0.0036) 

4)1,4)p, 4>(4) = (0.52,0.58,3.95,3.47)rad 
(toi6, Mi/ 2 , To,/x(Mz)) = (25000,280,-51380,1212) GeV 
{{mHi/miQY, (TOH„/mi6)^, tan/3) = (1.86,1.61, 50.29) 

= (9.2,578.8,35054.2) x lO^GeV 


Observable 

Fit 

Exp. 

Pull 

(7 

Mz 

91.1876 

91.1876 

0.0000 

0.4541 

Mw 

80.4784 

80.3850 

0.2320 

0.4027 

1 /<Teni 

137.2810 

0.0073 

0.3569 

0.6864 

X 10 ^ 

1.1789 

1.1664 

1.0598 

0.0118 

az{Mz) 

0.1192 

0.1185 

0.8199 

0.0008 

Mt 

174.0947 

173.2100 

0.7171 

1.2337 

mb{mb) 

4.1986 

4.1800 

0.5092 

0.0366 

m-T 

1.7772 

1.7768 

0.0428 

0.0089 

Mb-M, 

3.1701 

3.4500 

0.8720 

0.3209 

mc(mc) 

1.2509 

1.2750 

0.9333 

0.0258 

m,(2 GeV) 

0.0953 

0.0950 

0.0609 

0.0050 

md/m ,(2 GeV) 

0.0702 

0.0513 

2.8247 

0.0067 

to 

0.0018 

0.0019 

0.4528 

0.0001 

Mp 

0.1056 

0.1057 

0.0578 

0.0005 

Me X 10'‘ 

5.1143 

5.1100 

0.1674 

0.0256 

IKel 

0.2245 

0.2253 

0.5931 

0.0014 

Ted 

0.0404 

0.0408 

0.1670 

0.0021 

IKbl X 10^ 

3.1235 

3.8500 

0.8446 

0.8601 

TmI X 10=* 

8.8463 

8.4000 

0.7418 

0.6016 

ITd 

0.0396 

0.0400 

0.1508 

0.0027 

sin 2 P 

0.6296 

0.6820 

2.7214 

0.0193 

I^K 

0.0022 

0.0022 

0.0022 

0.0002 

AMbJAMb, 

34.8195 

35.0345 

0.0308 

6.9747 

AMbj X 10 *=* 

3.9946 

3.3370 

0.8224 

0.7996 

X 10=> 

7.5883 

7.5550 

0.0621 

0.5363 

mgi X 10^ 

2.4649 

2.4620 

0.0197 

0.1455 

sin=* 012 

0.3028 

0.3070 

0.1125 

0.0370 

sin=* 023 

0.6600 

0.5125 

1.1300 

0.1305 

sin=* 013 

0.0162 

0.0218 

1.7510 

0.0032 

Mh 

126.2697 

125.7000 

0.1882 

3.0265 

BR(B 57 ) X 10^ 

2.7220 

3.4300 

0.5419 

1.3064 

BR(Bs — 1- X 10^ 

2.7213 

2.8000 

0.0888 

0.8867 

BR{Bd — 7- X 10^*^ 

1.0734 

3.9000 

1.7509 

1.6143 

BR(B — f Ti/) X 10^ 

6.2223 

11.4000 

1.3588 

3.8104 

BR{B —> )i< 52<6 GeV^ ^ 10^=* 

4.7860 

3.4000 

0.2739 

5.0610 

BR(B —> K*^ )i 4.18<92<16 GeV 2 X 10® 

7.5495 

5.6000 

0.1356 

14.3788 

0 o(Mb(B ^ K*^l+tl-)) 

3.7120 

4.9000 

0.9190 

1.2927 

PhiB -> K* GeV 2 

0.7207 

0.6500 

0.2101 

0.3366 

Fl{B K*fX )i4.i8<52<i 6 GeV 2 

0.3108 

0.3300 

0.0726 

0.2644 

P2(B -4 ii'*/4+/4“)i<,j2<g GeV^ 

0.0331 

0.3300 

2.3939 

0.1240 

P2(B -4 K*I2'*'I2 )i4.i8<g2<i6 GeV2 

-0.4336 

-0.5000 

0.3364 

0.1974 

Pl{B -4 ii'*/4'''/4“)i<.j2<g GeV= 

0.5717 

0.5800 

0.0208 

0.3988 

Pl{B ^ K*fl )i4.i 8<?2<16 GeV 2 

1.2190 

-0.1800 

1.7066 

0.8198 

Pi(B “4 K*11 )l<,j 2 <e GeV 2 

-0.4335 

0.2100 

2.2451 

0.2866 

Pi(B —> K*fX )i4.i8<52<i 6 GeV 2 

-0.7116 

-0.7900 

0.1552 

0.5052 

Total 

46.7692 
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Table 7: Benchmark point with mie = 25 TeV, = 1.187 TeV, a = 1.5: 
(l/acMcea) = (25.95,2.60 x lO^® GeV,-1.50%) 

(A, Ae, a, Ae, p, Ae', Ae^ = (0.6100,0.0310,1.1459,0.0049, 0.0708, -0.0019,0.0037) 
((/)^, (/>,-, (l)p, (j)^) = (0.53, 0.57, 3.94, 3.49)rad 
(mi6, Mi/ 2 , To, M(Tfz)) = (25000,520,-51157,1236) GeV 
{{mHa/rnie)‘^,{mH^/rnie)‘^,tan/3) = (1.85,1.61,50.19) 

{Mr^,Mr^,Mr^) = (9.1,567.0,32370.5) x lO^GeV 


Observable 

Fit 

Exp. 

Pull 

(7 

Mz 

91.1876 

91.1876 

0.0000 

0.4540 

My\r 

80.5197 

80.3850 

0.3344 

0.4029 

l/c^ein 

137.1416 

0.0073 

0.1540 

0.6857 

X 10^ 

1.1829 

1.1664 

1.3978 

0.0118 

a^(Mz) 

0.1189 

0.1185 

0.4798 

0.0008 

Mt 

173.8449 

173.2100 

0.5150 

1.2328 

mb{mt) 

4.2023 

4.1800 

0.6094 

0.0366 

TTlr 

1.7772 

1.7768 

0.0450 

0.0089 

Mb-Me 

3.1680 

3.4500 

0.8791 

0.3207 

mc{mc) 

1.2570 

1.2750 

0.6979 

0.0258 

ms{2 GeV) 

0.0947 

0.0950 

0.0671 

0.0050 

md/ms{2 GeV) 

0.0700 

0.0513 

2.7901 

0.0067 

1/Q' 

0.0018 

0.0019 

0.5027 

0.0001 

Mp 

0.1056 

0.1057 

0.1457 

0.0005 

Me X lO'* 

5.1145 

5.1100 

0.1775 

0.0256 

IKsI 

0.2244 

0.2253 

0.6440 

0.0014 

IVfcl 

0.0407 

0.0408 

0.0584 

0.0021 

\Vnb\ X 10^ 

3.1307 

3.8500 

0.8363 

0.8601 

\Vtd\ X 103 

8.8596 

8.4000 

0.7639 

0.6016 

IV.I 

0.0398 

0.0400 

0.0652 

0.0027 

sin 2/3 

0.6285 

0.6820 

2.7790 

0.0193 


0.0023 

0.0022 

0.1149 

0.0002 

AMbJAMb, 

35.5946 

35.0345 

0.0786 

7.1295 

AMb^ X 10^^ 

3.9756 

3.3370 

0.8025 

0.7958 

rriji X 10^ 

7.6111 

7.5550 

0.1046 

0.5364 

rrigj X 10^ 

2.4657 

2.4620 

0.0255 

0.1455 

sin^ 012 

0.3134 

0.3070 

0.1724 

0.0370 

sin^ 023 

0.6319 

0.5125 

0.9146 

0.1305 

sin^ 013 

0.0153 

0.0218 

2.0337 

0.0032 

Mb 

124.5455 

125.7000 

0.3814 

3.0265 

BR{B — ^ sj) X lO'^ 

2.7270 

3.4300 

0.5372 

1.3087 

BR{Ba — ^ ) X 10® 

2.5215 

2.8000 

0.3228 

0.8627 

BR'{Bd Mm^) X 10“ 

1.0192 

3.9000 

1.7861 

1.6129 

BR{B tv) X 10® 

6.2272 

11.4000 

1.3568 

3.8124 

BR{B —^ )i<g2<g Qgv2 X 10* 

4.8580 

3.4000 

0.2839 

5.1361 

BR{B —^ K*R^H )i4.18<ii2<16 GeV^ ^ 10* 

7.6648 

5.6000 

0.1415 

14.5975 

ql{ABB{B^K*^l+^Ji-)) 

3.7150 

4.9000 

0.9163 

1.2933 

Fl{B —>■ K*fl+H )i<^2<g GeV^ 

0.7208 

0.6500 

0.2103 

0.3366 

Fl{B —>■ K*qi )i4,i8<,2<ig GeV^ 

0.3108 

0.3300 

0.0726 

0.2644 

P2(S —)l<g2<g GeV^ 

0.0335 

0.3300 

2.3879 

0.1242 

P2{B^K*R^^ )l4,18<g2<ig GeV^ 

-0.4336 

-0.5000 

0.3364 

0.1974 

P'iiB )i<g2<g GeV^ 

0.5697 

0.5800 

0.0258 

0.3985 

PiiB K*^l+fJ, )i4,i8<g2<ig GeV^ 

1.2190 

-0.1800 

1.7066 

0.8198 

PiiB K*fJ,+ll )i<g2<g GeV^ 

-0.4334 

0.2100 

2.2450 

0.2866 

PiiB^K*fl+fJ, )i4,i8<o2<ig GeV^ 

-0.7117 

-0.7900 

0.1550 

0.5052 

Total 

47.7692 



17 




Table 8: Benchmark point with mie = 25 TeV, Mg = 1.613 TeV, a = 0: 
(l/acMces) = (26.22,2.32 x lO^® GeV,-0.65%) 

(A, Ae, a, Ae, p, Ae', Ae^ = (0.6096,0.0311,1.1398,0.0049, 0.0710, -0.0019,0.0038) 
((/)^, (/>,-, (l)p, (j)^) = (0.53, 0.56, 3.95, 3.49)rad 
(mi6, Mi/ 2 , To, M(Tfz)) = (25000,450,-51341,1226) GeV 
{{mHa/rnie)‘^,{mH^/rniG)‘^,tan/3) = (1.86,1.61,50.30) 

{Mr^,Mr^,Mr^) = (9.1,572.4,32277.4) x lO^GeV 


Observable 

Fit 

Exp. 

Pull 

(7 

Mz 

91.1876 

91.1876 

0.0000 

0.4535 

My\r 

80.4507 

80.3850 

0.1633 

0.4025 

l/c^ein 

137.7125 

0.0073 

0.9825 

0.6886 

X 10^ 

1.1732 

1.1664 

0.5798 

0.0117 

a^(Mz) 

0.1188 

0.1185 

0.4140 

0.0008 

Mt 

174.1882 

173.2100 

0.7927 

1.2340 

mb{mt) 

4.1954 

4.1800 

0.4220 

0.0366 

TTlr 

1.7781 

1.7768 

0.1417 

0.0089 

Mb-Me 

3.1568 

3.4500 

0.9175 

0.3196 

mc{mc) 

1.2595 

1.2750 

0.5993 

0.0258 

ms{2 GeV) 

0.0939 

0.0950 

0.2147 

0.0050 

md/ms{2 GeV) 

0.0701 

0.0513 

2.8052 

0.0067 

1/Q' 

0.0018 

0.0019 

0.5139 

0.0001 

Mg 

0.1056 

0.1057 

0.1818 

0.0005 

Me X lO'* 

5.1145 

5.1100 

0.1749 

0.0256 

IKsI 

0.2244 

0.2253 

0.6763 

0.0014 

IVfcl 

0.0404 

0.0408 

0.1729 

0.0021 

\Vnb\ X 10^ 

3.1033 

3.8500 

0.8681 

0.8601 

\Vtd\ X 103 

8.8101 

8.4000 

0.6817 

0.6016 

IV.I 

0.0396 

0.0400 

0.1531 

0.0027 

sin 2/3 

0.6270 

0.6820 

2.8562 

0.0193 


0.0022 

0.0022 

0.2052 

0.0002 

AMbJAMb, 

35.3739 

35.0345 

0.0479 

7.0854 

AMb^ X 10^^ 

3.9433 

3.3370 

0.7681 

0.7894 

rriji X 10^ 

7.6562 

7.5550 

0.1886 

0.5364 

rrigj X 10^ 

2.4631 

2.4620 

0.0077 

0.1455 

sin^ 012 

0.3170 

0.3070 

0.2689 

0.0370 

sin^ 023 

0.6264 

0.5125 

0.8722 

0.1305 

sin^ 013 

0.0149 

0.0218 

2.1658 

0.0032 

Mb 

124.5054 

125.7000 

0.3947 

3.0265 

BR{B — ^ sj) X lO'^ 

2.6840 

3.4300 

0.5789 

1.2887 

BR{Ba — ^ ) X 10® 

3.0247 

2.8000 

0.2429 

0.9252 

BR'{Bd Mm^) X 10“ 

1.1022 

3.9000 

1.7323 

1.6151 

BR{B tv) X 10® 

6.1884 

11.4000 

1.3727 

3.7966 

BR{B — ^ )i<«^<6 GeV^ ^ 10* 

4.7640 

3.4000 

0.2707 

5.0381 

BR{B — ^ K*M )i4.18<ii 2<16 GeV^ ^ 10* 

7.5110 

5.6000 

0.1336 

14.3059 

ql{ABB{B^K*^l+^^-)) 

3.6690 

4.9000 

0.9579 

1.2850 

Fl{B —>■ K*M )l<g2<6 GeV^ 

0.7225 

0.6500 

0.2149 

0.3374 

Fl{B —>■ K*M )i 4.18<«2<16 GeV^ 

0.3108 

0.3300 

0.0726 

0.2644 

P2{B K*MP‘ )l<g2<6 GeV^ 

0.0228 

0.3300 

2.5196 

0.1219 

P2{B^K*MfJ' )l4.18<g2<16 GeV^ 

-0.4336 

-0.5000 

0.3364 

0.1974 

Pi{B /f*/r+/i )i<g2<g GeV^ 

0.5820 

0.5800 

0.0050 

0.4001 

PiiB )l4.18<g2<16 GeV^ 

1.2190 

-0.1800 

1.7066 

0.8198 

PiiB K*MP‘ )l<g2<6 GeV^ 

-0.4455 

0.2100 

2.2578 

0.2903 

PiiB ^ K* fl+fJ- )l4,l8<o2<ie GeV^ 

-0.7116 

-0.7900 

0.1552 

0.5052 

Total 

48.8413 
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Table 9: Benchmark point with mie = 25 TeV, = 1.690 TeV, a = 1.5: 
(l/acMcea) = (26.38,2.09 x lO^® GeV,0.02%) 

(A, Ae, a, Ae, p, Ae', Ae^ = (0.6096,0.0311,1.1384,0.0049, 0.0708, -0.0019,0.0037) 
((/)^, (/>,-, (l)p, (j)^) = (0.52, 0.56, 3.96, 3.49)rad 
(mi6, Mi/ 2 , To, M(Tfz)) = (25000,900,-50846,1529) GeV 
{{mHa/rnie)‘^,{mH^/rniG)‘^,tan/3) = (1.86,1.60,50.31) 

{Mr^,Mr^,Mr^) = (9.1,579.0,32367.3) x lO^GeV 


Observable 

Fit 

Exp. 

Pull 

(7 

Mz 

91.1876 

91.1876 

0.0000 

0.4540 

My\r 

80.4655 

80.3850 

0.2000 

0.4026 

l/c^ein 

137.7323 

0.0073 

1.0111 

0.6887 

X 10^ 

1.1740 

1.1664 

0.6469 

0.0117 

a^(Mz) 

0.1188 

0.1185 

0.2979 

0.0008 

Mt 

174.3427 

173.2100 

0.9175 

1.2345 

mb{mt) 

4.2001 

4.1800 

0.5479 

0.0366 

TTlr 

1.7774 

1.7768 

0.0644 

0.0089 

Mb-Me 

3.1659 

3.4500 

0.8863 

0.3205 

mc{mc) 

1.2574 

1.2750 

0.6825 

0.0258 

ms{2 GeV) 

0.0936 

0.0950 

0.2741 

0.0050 

md/ms{2 GeV) 

0.0701 

0.0513 

2.8082 

0.0067 

1 /Q' 

0.0018 

0.0019 

0.5170 

0.0001 

Mp 

0.1056 

0.1057 

0.1571 

0.0005 

Me X lO'* 

5.1139 

5.1100 

0.1545 

0.0256 

IKsI 

0.2244 

0.2253 

0.6688 

0.0014 

Tcfcl 

0.0400 

0.0408 

0.3609 

0.0021 

\Vnb\ X 10^ 

3.0662 

3.8500 

0.9113 

0.8601 

\Vtd\ X 103 

8.7156 

8.4000 

0.5247 

0.6016 

IT.I 

0.0392 

0.0400 

0.2960 

0.0027 

sin 2/3 

0.6259 

0.6820 

2.9122 

0.0193 


0.0022 

0.0022 

0.0834 

0.0002 

AMbJAMb, 

34.7964 

35.0345 

0.0342 

6.9701 

AMb^ X 10^^ 

3.8958 

3.3370 

0.7165 

0.7799 

rriji X 10^ 

7.6614 

7.5550 

0.1984 

0.5364 

rrigj X 10^ 

2.4606 

2.4620 

0.0094 

0.1455 

sin^ 012 

0.3197 

0.3070 

0.3423 

0.0370 

sin^ 023 

0.6197 

0.5125 

0.8210 

0.1305 

sin^ 013 

0.0146 

0.0218 

2.2520 

0.0032 

Mb 

122.0502 

125.7000 

1.2059 

3.0265 

BR{B — ^ sj) X lO'^ 

2.6310 

3.4300 

0.6321 

1.2640 

BR{Ba — ^ ) X 10® 

3.5145 

2.8000 

0.7203 

0.9920 

BR'{Bd Mm^) X 10 “ 

1.0522 

3.9000 

1.7647 

1.6138 

BR{B tv) X 10® 

6.1009 

11.4000 

1.4090 

3.7610 

BR{B —^ )i<g2<g Qgv2 X 10* 

4.6780 

3.4000 

0.2583 

4.9484 

BR{B — ^ K*R^H )i4.18<ii2<16 GeV^ ^ 10* 

7.4066 

5.6000 

0.1281 

14.1080 

ql{ABB{B^K*^l+^Ji-)) 

3.6290 

4.9000 

0.9946 

1.2779 

Fl{B —>■ K*fl+H )i<^2<g GeV^ 

0.7240 

0.6500 

0.2189 

0.3380 

Fl{B —>■ K*qi )i4,i8<,2<ig GeV^ 

0.3108 

0.3300 

0.0726 

0.2644 

P2{B )l<g2<g GeV^ 

0.0132 

0.3300 

2.6254 

0.1207 

P2{B^K*R^^ )l4,18<g2<ig GeV^ 

-0.4337 

-0.5000 

0.3358 

0.1975 

P'iiB )i<g2<g GeV^ 

0.5918 

0.5800 

0.0294 

0.4014 

PiiB K*^l+fJ, )i4,i8<g2<ig GeV^ 

1.2190 

-0.1800 

1.7066 

0.8198 

PiiB K*fJ,+ll )i<g2<g GeV^ 

-0.4562 

0.2100 

2.2685 

0.2937 

PiiB^K*fl+fJ, )i4,i8<o2<ig GeV^ 

-0.7117 

-0.7900 

0.1550 

0.5052 

Total 

52.6056 
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B Natural SUSY Breaking 

Consider Heterotic orbifold models with dilaton/moduli SUSY breaking as discussed in [15] . 
Following Eqn. (60, 61) in this paper we have the scalar masses for sparticle a, the trilinear 
couplings and the gaugino masses are given by 

= ^ 3 / 2(1 + cos^ Ofia ■ 0^) + Vq (11) 

6 

Ao = -V3Cm3/2{sm6 *6*^^ + cos 0 ^ [1 + nl* + nf + - (T* + T*)di log Yigioie]), 

i=l 

( 12 ) 

Ma — Vlm 3/2 sin , (13) 

where 




1 + 


K) 

3m?o ’ 


(14) 


Vo is the tree-level cosmological constant that we will set to be zero, i.e. C = 1 and 6 
determines the amount of SUSY breaking of the dilaton sector versus the moduli sector, 
cos 6* = 0 means all the SUSY breaking is due to the dilaton. are the modular weights of 
matter helds and 0 gives the probability for the SUSY breaking contribution of each modulus, 
such that XlLi ~ then let 6^ = 0 such that at tree level gauginos are massless. 

They would then obtain one loop suppressed masses due to moduli SUSY breaking. We also 
assume that the Yukawa couplings, are independent of the moduli. As a particular 
example, consider a Z 2 ® Zg orbifold with twist vectors given by (1/2,1/2,0), (1/6,2/3,1/6) 
in the three two torii [16]. There are three Kahler moduli in this example. Then we have 


m 


2 

a 


= + 3(nje5 + n?e|)) 


(15) 


where we assumed that 02 = 0 . 

We now assume that the Higgs 10-plet comes from the bulk on the second two torus. 
Thus UiQ = = 0, = —1. Hence, mio = ^ 3/2 and therefore, mfg = ^^ 0(1 + 

If the 16-plet lives in the fifth twisted sector with modular weights, 
nfg = = —1/6, = —2/3, we have m^g = or mio = \/2miQ. This is our first 

constraint, that mio = (^m + ^jrd )/2 ~ A/ 2 mi 6 . 

We then have A(i 6 ig 10 ) = -AAmi6(e“*'^i0i[l-F 2 ? 7 ,fg+ e“*^ 303 [l-F2n^^and 
we need A(i 6 ig 10 ) = ^0 = — 2 mig where the last equality is the boundary condition at the 
sweet spot. This can be solved with e“*'>'i0i = e “*'^303 = and 7 = 7 r/ 6 . Of 

course, it would require the dynamics of stabilizing moduli and SUSY breaking to fix these 
particular values of 0 *. 
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